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Abstract. The Wbb, Zbb and ttbb production at LHC are irreducible backgrounds for possible observabil- 
ity of the Standard Model and Minimal Supersymmetric Standard Model light Higgs boson in processes 
involving associated WH, ZH and ttH production followed by the H —* bb decay. The comparison pre- 
sented in this paper uses the background estimates obtained with (a) complete massive matrix element 
implemented in AcerMC Monte Carlo generator and (b) PYTHIA implementation of the inclusive W , Z, and 
tt production, followed by the parton showering mechanism. Both approaches lead to a production of the 
final state of interest but differ in the approximations used. The focal point of this study is the comparison 
of the two approaches when estimating the background contributions to the light Higgs boson searches at 
the LHC. 



1 Introduction 

The Wbh, Zbb and ttbb production at LHC has been rec- 
ognized , see e.g. to present the most substantial ir- 
reducible backgrounds for the Standard Model (SM) and 
Minimal Supersymmetric Standard Model (MSSM) light 
Higgs boson observability in the associated production, 
namely WH, ZH and ttH, followed by the H ^ bb de- 
cay. The 'light Higgs boson' is in this context understood 
as the Higgs boson having the mass between 90 GeV and 
130 GeV, thus describing the SM and MSSM Higgs bo- 
son(s) in the mass range indicated hy the excess observed 
in searches of LEP experiments lH,^. 

The potential of the ATLAS detector at LHC for the 
SM and MSSM Higgs boson observability in the ttH pro- 
duction has already been carefully studied and documen- 
ted in PI and H . The proposed analysis requires four iden- 
tified (tagged) b- jets, reconstruction of both top-quarks 
decaying in hadronic and leptonic modes and an observ- 
able peak in the invariant mass distribution of the re- 
maining b-jet pair. The irreducible ttbb background is es- 
timated to contribute about 60-70% to the total back- 
ground, which consists mostly of the processes with a tt 
pair participating in the final state. The expected signifi- 
cance is deemed to be of the order of 3.6 a for the Higgs- 
boson mass of 120 GeV at integrated luminosity of 30/5~^, 
with the expected signal-to-background ratio in the mass 
window rated to be 32%. The total contribution from the 
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Wjjjjjj backgrounds was estimated to be an order of 
magnitude smaller than the one consisting of tijj events 
and the contribution of the ttZ production process was 
estimated to be negligible. 



The potential for the Higgs boson observability in the 
WH {W — > (i') production, considered as rather weak, 
is documented in ||l[ and ||^. The expected production 
rates would be sufficient for the signal observability in 
the mass range around 120 GeV only if the backgrounds 
from Wbb and ti events could be efficiently suppressed. 
For this channel both b-quarks are required to be tagged 
as b-jets and the reconstruction of a peak in the invari- 
ant mass of the b-jet system concentrated in the interval 
± 20 GeV around the expected Higgs mass could lead 
to an evidence of the signal. The identification of the ac- 
companying lepton is also required in order to trigger the 
data acquisition. Relatively simple topology of the final 
state does not leave much room for a severe optimisa- 
tion of kinematic cuts. The only possibility which can be 
explored seems to be using a veto on an additional jet 
to suppress reducible backgrounds or angular correlations 
between reconstructed b-jets themselves and/or between 
b-jets and leptons._ Within the low luminosity operation, 
the irreducible Wbb background contributes about 35% of 
the total background to this channel, the expected signal- 
to-background ratio in the given mass window being on 
the order of a few percent. The expected significance is 
estimated to be on the level of 3.0 a for the Higgs-boson 
mass of 120 GeV and integrated luminosity of 30/6"-'^. 



2 



B. P. Kersevan, E. Richter-Was: What is the Wbb, Zbb or ttbb irreducible background 



Much less promising, if not hopeless to consider, is 
the observability of the ZH production with subsequent 
leptonic Z-boson decay. Nevertheless, for the sake of com- 
pleteness, such a study was documented in Both b- 
quarks are required to be tagged as b-jets and the accom- 
panying leptons are required to be identified in order to 
trigger the data acquisition as well as for a reconstruction 
of the resonant peak around the Z-boson mass. The pres- 
ence of a peak in the reconstructed invariant mass of the 
b-jet system concentrated in a mass window of ± 20 GeV 
around the expected Higgs mass would lead to an evidence 
of the signal. The irreducible Zbb background contributes 
about 80% of the total background to this channel, the 
expected signal-to-background ratio in the mass window 
being on the order of a few percent. The expected signifi- 
cance is estimated to be about 1.0 cr for the Higgs-boson 
mass of 120 GeV and integrated luminosity of 30/6^^. 

It is thus evident that good theoretical understanding 
of the irreducible tibh, Wbb and Zbb backgrounds is crucial 
for the light Higgs boson observation at LHC. 

In the presented study the two available, albeit qual- 
itatively different, approaches for simulating these back- 
grounds are compared. The first one (ME) is to use the 
lowest order massive matrix elements of the AcerMC gen- 
erator which lead to the required final state. The latter 
is subsequently completed with the initial and final state 
radiation simulated via parton showering as implemented 
in the PYTHIA or HERWIG ||| generators. The second one 
(PS) is to simulate inclusive tt, W and Z production using 
the native PYTHIA or HERWIG implementations and subse- 
quently obtain accompanying b-quarks using the parton 
shower approximation only. Both approaches have their 
caveats. Using the complete 2 — > 4 matrix element might 
not be sufficient as the b-quarks could appear at the dif- 
ferent steps in the evolution of the partonic cascade and 
not necessarily only at the hard process level. Resorting 
to the parton shower approach on the other hand tends 
in several cases to underestimate hardness of the radiated 
partons and does not reproduce well their topological con- 
figurations. 

In the presented comparison the approach of concen- 
trating not on the partonic distributions but on the recon- 
structed experimental quantities, i.e. jets and isolated lep- 
tons, is chosen. The generated events are thus treated with 
a simplified reconstruction procedure using the algorithms 
of the fast simulation of the ATLAS detector at LHC 
and subsequent fiducial cuts on reconstructed jets and iso- 
lated leptons are applied roughly as foreseen for this type 
of physics at LHC detectors: geometrical acceptance for 
b-jets and isolated leptons identification down to pseudo- 
rapidity 77 of 2.5 and transverse momenta threshold for jets 
and leptons of 15 GeV. Although leptons will be identified 
with the lower thresholds, the 15 GeV transverse momenta 
represents roughly what is needed for triggering on such 
events. This limits the given comparison to the topologi- 
cal configurations close to those which will be selected by 
the experimental analysis. In jet reconstruction a simple 
procedure in form of a cone algorithm is performed, us- 
ing a cone radius of 0.4 in the pseudorapidity - azimuthal 



angle plane. Subsequently, the procedures of jet energy 
calibration and b-tagging are applied. In case jet-veto was 
also applied, events with additional jets having transverse 
momenta pr > 30 GeV and I77I < 5.0 were rejected. Ef- 
ficiencies for b-tagging and lepton identification are not 
included in the given numbers, only the efficiencies for jet 
and b-jet reconstruction are taken into account. More de- 
tails on t he p erformance of the applied algorithms can be 
found in |l|. 

The AcerMC Monte Carlo generator code and its in- 
terfaces to PYTHIA generator were used in the given eval- 
uation. The generated statistics was typically about 10^ 
events for the ME and 10* events for the PS simulation 
chain. The ME events were generated with AcerMC matrix 
element implementations and the PS events were gener- 
ated with the default settings of PYTHIA 6.2. The ME 
events were further completed with the initial and final 
state radiation to assure more realistic jet reconstruction 
efhciencies and multiplicities, thus leading to a better de- 
scription of the interaction kinematics. The CTEQ5L[pT| 
parton density functions were used for all estimations and 
proton-proton collision at 14 TeV centre-of-mass energy 
were simulated. A similar study could also be repeated 
using HERWIG instead of PYTHIA generator. 

At this point it might be relevant to mention that for 
the inclusive W and Z/j* production the so called im- 
proved parton shower algorithm is implemented in PYTHIA, 
i.e. some higher-order corrections are integrated; as shown 
in ||l^ , it gives good description of the complete pj^ spec- 
trum. A corresponding improvement is also expected for 
the production of jets in association with the W boson. 
For the tt production the standard parton algorithm is 
still used PYTHIA. 

The comparisons between these two generation ap- 
proaches for Wbb, Zbb and tibb events are discussed in 
Section 2, 3 and 4 of the paper. Final conclusions are sum- 
marised in Section 5. 



2 The Wbb irreducible background 

In this section the irreducible Wbb background to the 
Higgs boson searches in the WH production, followed by 
the H ^ bb decay, is discussed. The evaluation is based 
on two simulation approaches, ME and PS, as specified 
below: 

— ME: Use the 2^4 matrix element for qq — > M^(— > 
ii')g* {-^ bb) process as implemented in ^ . The a x BR 
= 33.2 pb for single leptonic flavour W-boson decay. 
The g* — > bb splitting is coded already at the level 
of the matrix element. This matrix element represents 
the lowest order contribution to the Ibb final state. 
The initial (ISR) and final state radiation (FSR) are 
simulated with parton shower of PYTHIA, followed by 
hadronisation in order to complete the event genera- 
tion. 

— PS: Use the 2 — > 1 matrix element for qq -^W process 
as implemented in PYTHIA, followed by the ISR/FSR 
and hadronisation. The a x BR = 17200 pb for the 
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Fig. 1. The transverse momenta of the W-boson. Solid line de- 
notes the PS events, dashed the ME events. Events were filtered 
as specified in the figure. 



single flavour leptonic W-boson decay. Gluon splitting 
in the ISR partonic cascade is the source of b-quarks in 
the event. The implementation of this process includes 
ISR/FSR modeled with an improved parton shower ap- 
proach |12,|l^ to match/merge with higher order ma- 
trix element calculations for the spectra. 



Considering the heavy flavour content of the cascade, 
a part of the higher order corrections is a priori already in- 
cluded in the parton shower approach (c.f. |T4|p^ ) whereas 
only the lowest order term for the g* bb splitting is 
present in the matrix element calculations. 

However, as the 2^4 matrix element with activated 
ISR and FSR is used, a part of the higher-order corrections 
is also to some extend included (e.g. additional branchings 

01 b bg are made possible), however without a rigorous 
check on the possible double- counting [p^ . In order to limit 
the double counting possibility the scale of the ISR/FSR 
shower is tuned to avoid the radiated partons being harder 
than the hard-process ones. What is not included in the 

2 — > 4 matrix element calculations is the contribution from 
events where gluon and quark interact in the hard process 
to produce the VF-boson, with the gluon splitting into b- 
quarks occurring in the further steps of the cascade. As 
well not included is the contribution where the final state 
gluon in the hard process decays into light quarks, while 
b-quarks appear in another branch of the partonic shower. 

Comparison of the differential distributions of the pjf 
spectra is presented in Fig. |l|. In the top plot the PS dis- 
tribution after filtering on events with reconstructed 2 jets 
-I- 1 lepton is shown, along with PS and ME distributions 
after filtering on exactly 2 b-jets -I- 1 lepton. In the lower 



Table 1. Cross-section for the qq Ivbb ME and qq W 
production PS with W ^ Iv decay (single flavour). Efficien- 
cies for h-tagging and lepton identification are not included, 
only included are the efficiencies for jet reconstruction and b- 
jet tagging. 



Selection 


qq —* Ivbb 
ME 


qq -> W(-^ £iy) 
PS 


Generated a x BR 


33.2 pb 


17200 pb 


Two b-jets -I- one lepton 
mbb-jcts = 100 - 140 GeV 


1.46 pb 
0.16 pb 


3.10 pb 
0.23 pb 


Two b-jets -I- one lepton + 
jet- veto 

mtb-jctB = 100 - 140 GeV 


1.13 pb 
0.12 pb 


1.55 pb 
0.12 pb 



range of the spectra the PS events have different slope 
when requiring 2 jets or 2 b-jets; requiring 2 b-jets strongly 
suppresses the selection of events in the lowp^ range. Af- 
ter filtering on 2 b-jets the slope and normalisation of the 
PS and ME events agree relatively well, the normalisation 
ratio being on the level of 1.5-2.5. A closer look (bottom 
plots) indicates a substantial enhancement of PS events in 
the range p^ = 40 — 120 GeV. For much higher pj^ the 
ME predictions start to exceed the PS ones. 

Table || quantifies the expected cross-sections for the 
inclusive production after requiring a reconstructed £bb 
final state. The PS predictions turn out to be 50-100% 
higher than the ME ones but are still quite compatible 
for events with the invariant mass of the b-jet system in 
the range of interest. After requiring jet-veto, important 
for selection of the WH channel to suppress the tt back- 
ground, the ME predictions agree with the PS ones in the 
mass range of interest. The numbers in this table illus- 
trate that the slope of the invariant mass distribution of 
the b-jet system and jet multiplicities is quite different for 
different simulation approaches. 

For the sake of the evaluation consistency, factorisa- 



tion/renormalisation energy scale — mfy was used 
throughout the ME and PS event generation. The factori- 
sation/renormalisation scale dependence for this process 
is of the order of 20% at most, which was estimated by 
using a range of different definitions of the energy scale 
implemented in 

The fact that the overall normalisations differ, one can 
naively interpret as an effect indicating that the effective 
cascade branching into heavy flavour quarks is more in- 
tense in PS events, quantitatively by almost a factor of 
1.5-2.5 higher. The feature was discussed in more detail 
in the conclusions being that it rather reflects a sum 
of several effects, not just a simple enhancement in the 
effective cascade branching. 

In Fig. ^ the distributions relevant for the experimen- 
tal analyses are plotted for events with two reconstructed 
b-jets and one tagged lepton. One can observe that there 
is a signiflcant difference in the predicted cone separation 
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Fig. 2. Top: the cone separation between b-jets, transverse mo- 
menta of the individual b-jets; Bottom: transverse momenta of 
the b-jet system and the invariant mass distribution of the b-jet 
system. Solid line denotes the PS events, the dashed one ME 
events. 




(GeV) (GeV) 



Fig. 3. The invariant mass distribution after gradually adding 
selection requirements is shown. Solid line denotes the PS 
events, the dashed one the ME events. 

between the b-jets, the i?bb-jotE|5 The parton shower ap- 
proach predicts more events with a small cone separation. 
These events can be rejected if a threshold on that sep- 
aration is required. In fact, they are not contributing to 
the higher range of the invariant mass of the b-jet pair. 
As expected, the shapes of the transverse momenta of the 
individual b-jets and of the b-jet system are also harder 
for the ME events. Quite different is the slope of the in- 
variant mass distribution of the b-jet system; for PS events 

^ Cone separation is calculated as the separation in the plane 
of pseudorapidity (7;) and azimuthal angle (0), the -Rbb-jcts ~ 
^ (Z\(;ibb-jots)^ + (^^rybb-jots)^. 



the distribution falls down more rapidly for higher masses. 
In Fig. 1^ the invariant mass distribution after gradually 
adding selection requirements is drawn. A better agree- 
ment is observed when a large cone separation between 
b-jets and the jet-veto on the additional jets are required. 
This could be explained by the fact that after these re- 
quirements the topology of PS events comes closer to the 
ME one. As the ME events exhibit lower average multi- 
plicity for reconstructed jets than the PS events, the latter 
are suppressed stronger by the jet veto requirement and 
thus give cross-section predictions which are even lower 
than the ones for ME events. The invariant mass distri- 
bution for ME events however remains harder than the 
one for PS events. Please also note the effect of imposing 
-Rbb-jots > 0.7 on the shape of the invariant mass dis- 
tribution of the b-jet system, which shifts the distribution 
maximum to about 40-60 GeV (compare Fig. I and I). 

More detailed discussions on the quantitative differ- 
ences between PS and ME approaches for the qq Wiy-^ 
lv)g{-^ bb) process can be found in 

3 The Zbb irreducible background 

In this section the estimates for the irreducible Zbb back- 
ground to the Higgs boson searches in the ZH production, 
followed by the H bb decay, are discussed. The evalua- 
tion is based on two simulation approaches, ME and PS, 
as specified below: 

— ME: Use the 2 — > 4 matrix element for gg, qq 
Z/j*(^ ££)bb process as implemented in AcerMC. The 
a X BB = 26.2 pb for single leptonic flavour Z/j* de- 
cay and for the invariant mass of the lepton pair above 
60 GeV. This matrix element represents lowest order 
contribution to the Mbb final state. The initial (ISR) 
and final state radiation (FSR) is simulated with par- 
ton shower of PYTHIA, followed by hadronisation to 
complete the event generation. 

— PS: Use the 2^1 matrix element for qq Zj^* 
process as implemented in PYTHIA, followed by the 
ISR. The a x BR = 1640 pb for single lepton flavour 
decay and invariant mass of the lepton pair above 
60 GeV. Gluon splitting in the ISR partonic cascade 
is the source of b-quarks in the event. The implemen- 
tation of this process includes ISR/FSR modeled with 
an improved parton shower approach [p^ , ||l3|] which 
integrates some higher-order corrections. 

The study begins by comparing differential distribu- 
tions of the Z-boson transverse momenta for both simula- 
tion approaches, see Fig. ^, after requiring reconstructed 
leptons and jets (b-jets) in the final state. With the gener- 
ation threshold for the invariant mass of the lepton pair at 
60 GeV, the studied events are dominated by the on-shell 
Z-boson exchange. 

Contrary to the Wbb case, the slopes of the p^ distribu- 
tions in the ME and PS events are very different; the ME 
events are found to be much harder. One should remember 
that there are several topologies of the Feynman diagrams 
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Fig. 4. The transverse momenta of the Z-boson. Solid line 
denotes the PS events, dashed one the ME events. Events were 
filtered as specified in the figure. 

Table 2. Cross-section for the gg, qq Ubb ME and qq — > 
Z/^* PS events with Z/-y* II decay (single flavour). 
Throughout the study the mass of the lepton-pair is required 
to be above 60 Gc V. Efficiencies for b-tagging and lepton iden- 
tification are not included, only the efficiencies for jet recon- 
struction and b-jet tagging are taken into account. 



Selection 


9g,<ig-^ Ubb 
ME 


qq^Zh*{-^U) 

PS 


Generated: a x BR 


26.2 pb 


1640 pb 


Two b-jets + 
two leptons 


1.70 pb 


1.62 pb 


m« = ± 10 GeV 


1.54 pb 


1.48 pb 


n-bb-jcts = 

100 - 140 GeV 


0.28 pb 


0.31 pb 



leading to the Ubb final state, see e.g. the dominant 
one being the contribution from the multipheral topolo- 
gies of gg ^ Z/j*bb where each gluon splits into a bb pair 
and the bb pair originating from different gluons annihi- 
lates to produce the Z-boson. This might explain why in 
this case the universal improved parton shower approxi- 
mation implemented in PYTHIA for qq Z/^* process is 
not working as well as in the previous case of the qq ^ W 
process and the £bb events. 

The overall normalisation prediction of both simula- 
tion streams_, see Table seems to be in reasonable agree- 
ment for ££bb events. This agreement remains after requir- 
ing lepton-pair and b-jet pair within the mass window, but 



Fig. 5. Distributions of events with two tagged b-jets and a lep- 
ton pair with invariant mass around the mass of the Z-boson: 
The cone separation between b-jets, transverse momenta of in- 
dividual b-jets, and transverse momenta of the b-jet system. 
Solid line denotes the ME events, the dashed one the PS events. 



would deteriorate if experimental analysis became sensi- 
tive to the hard tail of the distribution. In both genera- 
tion approaches the factorisation/renormalisation energy 
scale = m% used. Also in this case, the variation 
of the cross-section with the factorisation/renormalisation 
scales implemented in AcerMC is less than 20%. 

In Fig.|5|the distributions relevant for the experimental 
analyses are drawn. One can observe that, contrary to the 
previous case, there is no significant difference in the cone 
separation between b-jets (i?bb-jots) distribution. In both 
simulation approaches the dominant fraction of events has 
a b-jet pair with a large cone separation. This marks the 
topology of Ubb process as quite different from £bb events. 
Furthermore, the shapes of the transverse momenta of the 
individual b-jets and of the b-jet system are in this case 
also quite similar (compatible) for ME and PS events. 

Surprisingly similar, given the complexity of the topo- 
logies introduced by the Feynman diagrams, is the distri- 
bution of the invariant mass of the b-jet system in PS and 
ME events. This is illustrated in Fig. |[ To illuminate this 
further. Fig. |^ shows the separate ME contributions from 
the qq and gg events to the total invariant mass spectrum 
of the b-jet system. One sees a distinctly different shape 
of both components; as expected the qq component is very 
similar to the one of Wbb events. Nevertheless, the latter 
class of events contributes only on the level of 10% to the 
total. 
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Fig. 6. The distribution of the invariant mass of the b-jet sys- 
tem for events with two b-jets and lepton pair within the Z- 
boson mass window is shown. Solid line denotes the ME events, 
the dashed one the PS events. 
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Fig. 7. The distribution of the invariant mass of the b-jet sys- 
tem for events with two b-jets and lepton pair within the Z- 
boson mass window, generated with ME, is shown. The gg and 
qq components are shown separately. 

4 The ttbb irreducible background 

In this section the estimates for the irreducible tibb back- 
ground to the Higgs-boson searches in the tiH production, 
followed by the H ^ bb decay, are discussed. Two genera- 
tion approaches, ME and PS, which lead to the ttbb final 
state are considered. 

— ME: Use the 2 — > 4 matrix element for gg, qq — > ttbb 
processes as implemented in For the QCD compo- 
nent the ax BR = 2.7 pb, with leptonic decay (electron 
or muon) of one W-boson and hadronic decay of the 
second one, both W bosons being produced in the top 
decays. For the EW component, gg, {Z/W/j* —>-)bb, 
the a X BR = 0.26 pb. These matrix elements rep- 
resent the lowest order contribution to the ttbb final 
state. The interference between QCD and EW compo- 
nent is not available in the present implementation of 
. In the numerical evaluation the central value of the 
factorisation scale ||l7|] , Qgr n = (toj -|-to///2)^, with 
rriH = 120 GeV was used for the QCD component and 
the = m?^ was used for electroweak one. For this 
process, as already stressed in different choices of 
the factorisation scale could lead to the cross-section 
estimates differing even by factor of four. Event gener- 



ation is completed by ISR/FSR and hadronisation as 
modeled in PYTHIA. 
— PS: Use the 2 — > 2 matrix element for gg,qq tt pro- 
cess as implemented in PYTHIA, followed by the ISR. 
The a x BR = 189 pb for leptonic decay (electron or 
muon) of one W-boson and hadronic decay of the sec- 
ond one, both W bosons being produced in the top 
decays. Gluon splitting in the ISR/FSR partonic cas- 
cade is the source of additional b-quarks in the event. 
The default factorisation energy scale of PYTHIA 6 . 2 
is used. 

There are two classes of processes which lead to the ttbb 
final state, the QCD and EW ones; for the corresponding 
Feynman diagrams see The PS events, where the hard 
process is just the top-quark pair production, contribute 
only the QCD component. Consequently, only the QCD 
ME component should be directly compared with the PS 
one. 

As an inclusive control distribution the transverse mo- 
menta spectra of the top quarks were chosen. Fig. ^ shows 
that there is quite a good agreement between PS and ME 
predictions. With the factorisation energy scale used for 
evaluating the ME predictions the absolute normalisation 
agrees within 20% and the ratio of the PS and ME distri- 
butions is amazingly flat. For other choices of the factori- 
sation energy scale, the ratio would be quite different (see 
Table with the total cross-sections in Q ) , but subsequent 
checks confirm that the distributions remain very similar. 

In the proposed experimental analysis Q both top- 
quarks have to be reconstructed and the remaining two 
b-jets are than considered as possible candidates for the 
Higgs boson decay products. In the present study this se- 
lection procedure was replaced by the tight matching re- 
quirements of the b-jets and their partonic origin. These 
lead to a very clean separation between b-jets originating 
from the top-quarks and those which are not. In what fol- 
lows only the distributions of the b-jets which are not 
identified as originating from the top-quark decays are 
considered. 

In Table || the expected total cross-sections and the 
cross-sections after simplified event selection are given. 
The EW ME component is on the level of 10% of the 
total ME cross-section and on the level of 25% of all ME 
events accepted in the mass window. The fact that the res- 
onant electroweak background is not negligible makes the 
prospects for the observability more difficult, especially 
for the Higgs boson masses closer to the mass of the Z- 
boson. This mass range has been already excluded for the 
SM Higgs boson [||, but for the MSSM scenarios is stiU 
in the region of a possible discovery While resonant 
EW component was estimated some time ago as negligi- 
ble in [Qq the non-resonant EW component was to our 
knowledge not considered prior to this study. It is clear 
that the evaluation presented in should be now revised 
to include the EW background properly. 

^ It was the consequence of the implementation for the QQZ 
process in PYTHIA 5 . 7 generator used at that time, which was 
not evaluating the total cross-section correctly. 
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Fig. 8. The transverse momenta of the top-quarks. Solid line 
denotes QCD PS events, dashed line the QCD ME events. 
Events were filtered as indicated on the plots. 

Table 3. Cross-sections for the QCD gg,qq ttbb, and EW 
gg ttbb and QCD gg, qq tt production PS with one W- 
boson from top-quark decaying leptonically (electron or muon), 
and the other one hadronically. 



Selection 


gg,qq ttbb 
(QCD ME) 


gg ttbb 
(EW ME) 


99, qq tt 
( QCD PS) 


Generated: 
ax BR 


2.7 pb 


0.26 pb 


189 pb 


4 b-jets + 

1 lepton + 

2 jets 


0.123 pb 


0.014 pb 


0.145 pb 


TTlbb-jcts = 

100-140 GeV 


0.013 pb 


0.003 pb 


0.014 pb 



Fig. 9. For events with four b-jets, isolated lepton and at least 
two light jets, the distributions for b-jets not originating from 
top-quark decays are drawn. Top: The cone separation between 
b-jets, transverse momenta of individual b-jets; Bottom: trans- 
verse momenta of the b-jet system. Solid line denotes the events 
generated with QCD ME for ttbb production, the dashed one the 
ones with QCD PS. 




m,,-im (GeV) m„.j„ (GeV) 



In Fig. ^ the distributions relevant for the experimen- 
tal analyses, simulated with QCD ME and PS simula- 
tion approaches, are drawn. One can observe some dif- 
ferences within the expected cone separation between b- 
jets, i?bb-jcts- In both simulation approaches the dominant 
fraction of events has the b-jet pair with a small cone sep- 
aration. Quite different are the shapes of the transverse 
momenta of individual b-jets and of the b-jet system. In 
particular, the distribution of the transverse momenta of 
the b-jet system is much harder in the ME events than in 
the PS events. Nicely enough, the shape of the invariant 
mass distribution of the b-jet system is quite similar in the 
relevant mass range, see Fig. |lO|. Moreover, the normalisa- 
tions are in a satisfactory agreement in the relevant mass 



Fig. 10. For events with four b-jets, isolated lepton and at 
least two light jets the invariant mass of the b-jet system for 
b-jets not originating from top-quark decays is drawn. Solid line 
denotes the events generated with QCD ME for ttbb production, 
the dashed one the events generated by QCD PS. 

range, for lower masses the PS predictions are exceeding 
the ME ones. 

So far only the comparison between ME and PS pre- 
dictions for the QCD ttbb events was discussed. The EW 
component can at present only be simulated with the 
ME implemented in AcerMC. The implementation gives 
a possibility to estimate either only the resonant part, the 
<?5i99 ~* ttZ/j* production with Z/j* bb decay, or 
to explore the full EW contribution, namely the process 
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Fig. 11. Invariant mass of the b-quarks system (top plot) and 
b-jet system (bottom plots) in EW ME gg — > tiZ{^ bb) and 
EW ME gg ~+ {Z/Wf-y* —*)ttbb events. The distributions are 
plotted only for b-quarks (resp. b-jets) originating in the hard 
process. 

gg tt {Z/W/j* ~*)bbtt. An implementation of the 
qq ^ tt ^ {Z /W/^* —*)bbtt is still missing, but will very 
likely contribute no more than 10-20% of the total EW 
background (assuming the same ratio as for QCD qq and 
gg contributions). 

Fig. |ll| shows the respective invariant mass distribu- 
tions of the b-quark pair (left plot) and the b-jet pair (mid- 
dle and right plots) not originating from top-quark decays 
in reconstructed tibb events, as estimated with either full 
or with only resonant EW ME processes. The presence of 
the flat non-resonant component, which is quite substan- 
tial with respect to the resonant one, is rather evident. 
One can also clearly see how the shape of the EW back- 
ground is smeared when going from b-quark distribution 
to b-jet distribution. 

Finally the QCD and EW component of the ME sim- 
ulation chain are added together. Fig. ^ again shows the 
respective invariant mass distributions of the b-quark pair 
(top plot) and b-jet pair (bottom plots) in reconstructed 
tibb events. Also shown separately are the QCD gg — > tibb, 
QCD qq tibb and EW gg {Z/Wh* -^)tibb compo- 
nents. The EW component contributes around 20-25% of 
the QCD one in the mass range of interest. 

It is crucial to notice that the EW component produces 
a resonant structure above the non-resonant QCD-t-EW. 
This will make a discovery of MSSM light Higgs in the tibh 
channel, with the mass between 95-110 GeV, even more 
difficult than assumed in One should however keep 
in mind that the relative contribution of QCD and EW 



Fig. 12. Invariant mass of the b-quarks system (top plot) and 
b-jet system (bottom plots) in QCD ME and EW ME events. 
The distributions are plotted only for b-quarks (resp. b-jets) 
originated from the hard process. 

components strongly depends on the chosen factorisation 
scale, see 0, and will thus require very good theoreti- 
cal understanding before the signal observability could be 
claimed at LHC in this channel. 



5 Conclusions 

In this paper a quantitative comparison between matrix 
element (ME) and parton shower (PS) approaches for gen- 
erating key irreducible backgrounds to the light Higgs bo- 
son searches at LHC with multi-b-jet final states was pre- 
sented. Three processes were discussed: the Wbb, Zbb and 
tibb production. The AcerMC Monte Carlo generator was 
used for simulating ME events while PYTHIA 6 . 2 was used 
for simulating PS events. 

— For the Wbb background there is a reasonable agree- 
ment of spectra in M events. Expected rates are 
however almost 1.5-2.5 higher in PS simulation. This 
enhancement comes mostly from configurations where 
b-jets are very close. Requiring large cone separation 
between b-jets and vetoing events with additional jets 
brings the ME estimates to be in acceptable agreement 
for the relevant invariant mass range of the b-jet sys- 
tem. 

— For the Zbb background the PS approach is clearly un- 
derestimating the hardness of the tail of the spec- 
tra. This is due to the dominant contribution from the 
multipheral Feynman diagrams, not reproduced well 
even with an improved parton shower implementation 
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in PYTHIA. The distributions of the invariant mass of 
the b-jet system and their normahsations are in rea- 
sonable agreement. 
— There are two components of the tibb background, the 
QCD and EW ones. Only the QCD component can 
be simulated with the PS approach. For QCD PS and 
QCD ME events there is a relatively good agreement in 
the shape and normalisation of the invariant mass dis- 
tribution of the b-jet system. The transverse momenta 
distribution of that system is however much harder 
in ME events. The EW ME component is not negli- 
gible and leads to the resonant structure in the total 
QCD-I-EW background. This would make observabil- 
ity of the mass range around the Z-boson mass more 
difficult than hoped so far. 

It is quite evident that well understood theoretical pre- 
dictions for these processes will provide the key for estab- 
lishing the Higgs signal observability; it is however difficult 
to draw universal conclusions from the comparisons pre- 
sented above. Therefore, it is rather encouraging that the 
PS and LO ME predictions are not very far off. The 50% 
differences e.g. for Wbb events in the overall normalisa- 
tion are still within expected uncertainties for this type of 
background estimates. It would nevertheless be very im- 
portant to perform a similar comparison with NLO ME 
predictions. Such implementations are already becoming 
available for Wbb events @ and Zbb events jl^. Nev- 
ertheless, at the time of the studies presented here, they 
were still not in the form allowing for straightforward ap- 
plicability. 

Given several possible tunings of the parton shower 
model as presently available in PYTHIA, one could proba- 
bly easily improve further on the agreement between both 
approaches. The framework prepared in AcerMC generator 
could well be a nice tool for such a tedious task. It is how- 
ever not clear, that the parameters should be tuned in a 
way to make the PS predictions agree with the ME ones. 
Rather, an enhanced theoretical understanding for which 
applications the PS or ME predictions are more credible 
and why should be achieved first. We hope that the results 
presented here could contribute to such discussions. 

Indeed, most of the ME/PS comparison studies so far 
concentrated predominantly on soft gluon emissions and/or 
NLO effects, where the hard processes under study are 
very simple and/or unproblematic (e.g. Drell-Yan Z bo- 
son production, ti production etc. . . ), whereas this study 
was concentrated on the limiting case of processes involv- 
ing either a large set of diagrams (from ME viewpoint) 
or quite hard gluon emission (from PS viewpoint). The 
latter case might be even more problematic since the soft 
(NLO) effects have already been (to some extent) properly 
incorporated into the PS algorithms jl^,^ whereas the 
hard gluon radiation (at high px) is a. feature that might 
be supposed not to work really well under PS approach 
by virtue of CLAP equations ||2^ and has received rela- 
tively little study so far. An additional point to support 
this claim is that the fiducial kinematic cuts applied in 
the above studies, which were chosen to reflect the cuts 
that will be applied in the future experimental analyses. 



to a certain extent suppress the soft (NLO) effects (e.g. 
cuts on minimal px, minimal cone separation etc.). This 
makes the discrepancies in the hard process description 
the focal point of attention. 

The issue of establishing a consistent procedure for the 
appropriate Monte Carlo generation |^ is crucial as the 
complete evaluation of the expected background should 
include its estimates for both irreducible and reducible^ 
components |^ . Quantifying the discrepancies between the 
ME and PS simulation approaches is important because 
it indicates what could be the expected systematic bias 
on the evaluation of the reducible backgrounds. The fact 
that PS and ME predictions are not very far off in the dis- 
cussed cases of the irreducible backgrounds is also encour- 
aging for using PS approach for simulating their reducible 
components. 

Results from the presented studies thus give us a stron- 
ger confidence in using the PS approach for simulating the 
complete backgrounds. The ME approach alone would be 
not sufficient but is nevertheless very valuable for quanti- 
fying the uncertainties of the PS approach in case of the 
irreducible background components. The results however 
also indicate that (for obvious reasons) the PS approach 
alone is also not without its shortcomings. Some contri- 
butions, like e.g. the full EW ttbb can presently not be 
covered by the PS algorithms. One should thus definitely 
aim for having both approaches available in a form which 
is straightforward to use in the experimental analyses with 
a clear theoretical understanding of their individual short- 
comings. 

This work was performed within the framework of the Higgs 
Working Group of the ATLAS Collaboration. We are grateful 
to all our colleagues for the inspiring atmosphere and several 
very valuable discussions. We have used a simplified version of 
the fast simulation of the ATLAS detector for the quantitative 
evaluations of the detector responses presented in this paper. 
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